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I. INTRODUCTION 

The studies of the cosmological aspects of supergravity 
and string theory have a long history, going back to the 
beginning of the 80's. At present, there is a new wave of 
interest in the cosmological aspects of string theory. The 
subject is rather complicated, partly due to the complex- 
ity of the analytical study of these theories. For example, 
to find the F-term part of the effective potential in super- 
gravity, one should specify the expressions for the Kahler 
potential K and superpotential W . Even with the sim- 
plest K and W, the computation of potential is tedious, 
especially if there is more than one superfield, and the 
resulting expression is hard to analyse. 

Deriving cosmological consequences from the models 
based on string theory and supergravity is intricate too. 
For non-canonical Kahler potentials (which are the rule 
rather than the exception) the equations of motion ac- 
quire additional velocity-dependent terms whose effects 
are not easily understood using intuition based on the 
simplest scalar field models. For the case of one field, 
one can always reduce the theory to the canonical form, 
but this method does not work for the description of 
the simultaneous motion of several different fields, so one 
should really solve the system of equations keeping the 
non-canonical kinetic terms throughout. 

Our Mathematica-based package "SuperCosmology" is 
intended to simplify the study of supergravity potentials 
and of the cosmological models based on string theory 
and supergravity. The programs we describe here were 
used in papers where only the final results of com- 

putations were presented. The purpose of this paper is 
to present the explanation of the programs and methods 
we used in computations. Also we present some new cos- 
mological models and use them to demonstrate how our 
package works. 



The SuperCosmology package consists of two parts. 
Part I has the following Mathematica nb-files: 

SuperPotential .nb, 
SuperPotential_KKLT .nb, 
SuperPotential J: ine_tune .nb, 
SuperPotential_D3 . nb. 

In SuperPotential . nb, SuperPotential_KKLT.nb, 
SuperPotential J: ine_tune .nb one finds examples from 
0, H and [jj, respectively, of using our program "Su- 
perPotential" in computation of a scalar potential. In 
SuperPotential_D3 .nb, a new example of an inflation- 
ary potential with a mobile D3 brane is studied, which 
is based on the D3/D7 inflationary model investigated 
before in [|-[l3|. 

In part II, we present the program "FRW" used in for 
the numerical solution of the Friedmann equations for a 
system with any number of scalar fields with geometric 
kinetic terms of the form ^Gij{(j), 4>*)d(j) l d(j)^ specified by 
a metric Gij{(j>) on the scalar manifold. Part II has the 
following Mathematica nb-files: 

FRW_N2.nb, 

FRW_DarkE.nb, 

FRW_LateDarkE.nb. 

In FRW_N2 .nb, we show a dark en ergy model based on 
the N=2 supergravity model 01 i uM which was also 
discussed in 

New results on the dark energy model based on the 
KKLT model are presented in FRWJDarkE.nb and 
FRW_LateDarkE.nb. The examples in FRW_DarkE.nb de- 
scribe the situation when the system has not yet reached 
the dS minimum (so that the scalars are still moving and 
flu is still increasing). In FRW_LateDarkE.nb we study 
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the long term evolution including the time when the dS 
minimum is reached by the scalars. 



II. N=l SUPERGRAVITY POTENTIALS 
FIXING THE MODULI 

It is believed that an N=l d=4 supergravity may serve 
as an effective theory for a more fundamental string/M 
theory. A generic supersymmetric gravity has two types 
of geometries: the space-time geometry 

ds 2 = g llv (x)dx fi dx u , 

and the moduli space geometry 

ds 2 = G l j(fadfadfa , 

which both affect the cosmological models. Here we con- 
sider a class of models of N=l supergravity with any 
number n c of chiral multiplcts (fa , fa — (fa)* , i — 
l,...,n c ), one Abelian vector multiplet and a gravita- 
tional multiplet. In application to cosmology, one pri- 
marily needs to know the total potential of the scalar 
fields and their kinetic terms (since these are typically 
non-canonical) . In the general case of many vector mul- 
tiplets, the kinetic function is a matrix f a p{fa) where 
a, f3 = 1, ...,n v with n v the number of vector multiplcts. 
Here we will consider a simpler picture with only one 
vector multiplet and one holomorphic function f(fa). 
The supergravity action is defined by the functions 0] 

W(4>), K(fafa), and f(fa). (1) 

(We do not include constant FI terms here, however, 
we will have a field-dependent D-term potential.) The 
bosonic part of the action is (we use units in which 
M P = 1): 

- i(Re/)F^^ + ii(Im/)F^F^. (2) 
The potential consists of an i^-term and a D-term: 

V = V F + V D , (3) 

V F = e K [(D t W){g- 1 ) l j (D 3 W*)-3WW*] , (4) 

V D = i (Re f)DD\ bos = ^(Ref^VV, (5) 

with 

V l W = d l W + {d l K)W , (6) 

V(fafa) = i [Sfad'Kifafa))] = 

= i [SpdiK^, <j>*))] , (7) 



where (5 fa, 6 fa) defines the £7(1) gauge transformations 
of chiral superfields. The covariant derivative of fa is 

d^fa = d^fa - W^Sfa . (8) 

Our package "SuperCosmology" is intended to help 
with studying models in the class presented above. In 
particular, it allows one to find the potential ©-© 
for a given input Q plus the gauge transformations of 
the superfields 8 fa. We will explain now how to use 
part I of our package by considering the example given 
in SuperPotential .nb. 

In the file SuperPotential .nb we use the scalar po- 
tential of the D3/D7 inflationary model with a light D7 
brane moving towards a heavy stack of D3 branes in an 
internal space with a stabilized volume 0. We give as 
input a list of complex fields p, S, <&+ , 3>_ . These are com- 
ponents of the holomorphic field fa. (In the files, we use 
barred symbols for the anti-holomorphic fields, e. g. p is 
the complex conjugate field to p.) The real parts of the 
fields p and S correspond to the volume of the internal 
space and the inflaton, respectively. The next piece of in- 
put is the holomorphic superpotential W(<j>), the Kahlcr 
potential K(fafa), and the U(l) gauge transformations 
of the scalar fields. The program calculates all impor- 
tant structures, including the cr-model metric Gi 3 (fa fa) 
defining the kinetic terms of the chiral superfields. It also 
presents a matrix form of it and an expression for all co- 
variant derivatives of the superpotential Di W. This sim- 
plifies the search for the supersymmetric configurations 
with constant scalar fields, since we can now easily solve 
the equation Di W = in all directions fa. Next, the pro- 
gram calculates the F-term and the D-term potentials as 
functions of all fields (fa fa). We found it convenient for 
computations to switch to real scalar fields by splitting 
the complex fields into their real and imaginary parts, 
p = a + ia, so that the total potential can be given as 
a function of 2n c real scalar fields. Finally, the program 
plots various sections of the potential as functions of just 
two fields at fixed values of other fields. The physical 
properties of this model are described in • 

We supplement SuperPotential .nb with a few more 
examples. The example in SuperPotentialJCKLT.nb 
shows how the stabilization of the volume was achieved 
in The example in SuperPotentialj£ine_tune .nb 
presents a detailed calculation of the fine-tuning proce- 
dure for the potential described in the Appendix F of 

Our fourth example presented in 
SuperPotential_D3 .nb, is based on the potential 
for the model of inflation with D3/D7 branes Q with 
volume stabilization, where the D7 brane is heavy and 
therefore only the D3 brane position modulus is turned 
on. Note that the Kahler potential in the case of a heavy 
stack of D3 branes and a light D7 0, which is used in 
SuperPotential .nb, is given by 

K=-3Hp + p)- iS ~ 2 S)2 . (9) 
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The opposite limiting case with a light D3 and heavy 
D7 branes, considered in SuperPotential_D3.nb, is de- 
scribed by a Kahler potential 



K = -3 In p 



(10) 



in agreement with 0, and the superpotential 

the same as in the KKLT model 



W = W + Ae~ ap . 



(11) 



The example given in the file SuperPotential_D3.nb 

describes a situation which was not explored before 
and turned out to be quite attractive. We have the 
same "trench" in the s — ^tdll direction which is to 
be interpreted as the inflaton direction. As shown in 



_ (p-p) 
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SuperPotential_D3 . nb, the dependence on a = 
quite simple: the potential includes ccosa with c neg- 
ative, and, therefore, there is a minimum at a = 0: 
the same behavior as in 0. However, the dependence 



on the field f3 



_ (4>+<t>) 



in this model is quite different 



from that in ■ Indeed, in the model studied in 
(see SuperPotential .nb), the F-term potential has a 
supersymmetric extremum at j3 = which is a saddle 
point. This saddle point remains after lifting the po- 
tential by adding a D-term, and a stable configuration 
must have some j3 ^ after (3 rolls down to the clos- 
est nearby minimum. In the model presented here in 
SuperPotential _D3 . nb, the F-term potential also has a 
saddle point at j3 = 0. However, here after the D-term 
potential is added, this saddle point turns into a perfect 
dS minimum. Since the stabilized value of the volume 
modulus a cr depends on the values of a and (3 at the ex- 
tremum, a cr here is the same as in the KKLT model 0, 
with a cr = (3 cr = 0. This gives a good starting point for 
a study of stringy models of hybrid inflation (in particu- 
lar, type IIB string theory compactified on K3 x ^ with 
D3/D7 branes present), based on the same mechanism of 
the volume stabilization as in the model of 41. 



III. SOLVING FRIEDMANN EQUATIONS IN 
SUPERGRAVITY-BASED MODELS OF DARK 
ENERGY 

In a series of papers 0-0, the issue of dark energy 
has been studied in the context of supergravity. In par- 
ticular, in these models we have studied the future of the 
universe and compared the predictions of some models 
with observations. The more recent supernova data for 
Z > 1 in £3 were also compared with the theoretical 
predictions in 0, 0- 

The action of a generic four-dimensional gauged su- 
pergravity which can be used for a dark energy hidden 
sector, includes gravity coupled to 2n scalar fields (f> 1 , and 
a potential: 

g-V 2 L = ~R + I GiM d»ftd v ft <T - V{ft) • (12) 



Here Gy {cf>) is the metric on the moduli space, related to 
the Kahler potential. (In the cosmological context it is 
convenient to work with 2n real scalar fields, as in l|12(l . 
rather then n complex fields). 

We consider here a model with dark energy represented 
by the scalars 4> z with the Lagrangian (|12J) . Also, we 
include the usual cold dark matter in the model, with 
the energy density pm, 



C 

Pm — q 
a 6 



(13) 



Also, we assume that the space is a flat FRW universe, 
ds 2 = dt 2 — a(t) 2 dx 2 , and that the scalar fields are homo- 
geneous. With these assumptions, the scalar and Fried- 
mann equations are: 

ft + 3-ft + T) k &^ k + G ij 7777 = , (14) 



a _ -pm + 2V ~ 4E ki . 
a 6 



(15) 



Here Tj k (<j)) are the Christoffel symbols in the moduli 
space defined by the metric (</>), and E k i n is the ki- 
netic part of the dark energy, 

1 



E, 



kin 



:G, 



1 M 



(16) 



A. First order Friedmann equations 

The Friedmann equations for homogeneous scalar 
fields take a particularly simple form in the Hamiltonian 
(first order) formulation, when the metric is a spatially 
flat FRW. For numerical calculations the first order equa- 
tions are also much more suitable so we use them in our 
computations. 

With the Lagrangian 11211 the canonical momenta are 



Pi =a 3 (t)G ij {tf>)<P' 1 , 

and the Hamiltonian is 
1 



(17) 



H(P,<j>,t) 



2a 3 (t) 



GPiflPiPj+artyVW) . (18) 



The equations of motion have the canonical form, 

■ i= m p __dH 
" op ' l ~ ■ 

Using the Hamiltonian l|18|) and adding the equation for 
the scale factor, we get a system of the first order ODEs: 



(19) 



(20) 



P 



1 dG 



kl 



2a 3 dft 
a = a H 



P k Pi - a' 



;dV(<t>) 



H 



-pM + 2U — 4:E k in 2 
tl 



(21) 
(22) 

(23) 
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Here the kinetic energy Ekin is given by 



1 



1 



E kin = 1 ^G«{4>)P i P j = ~G ij 



(24) 



and pm is given in |Q . 

In the simplest case of one dark energy scalar field with 
a minimal kinetic term of the form (4>) 2 /2 the system of 
equations pU)l. (|2*T)l reduces to a familiar form, 



-P 



P = -a 3 V' 



(25) 



< + 3-<j) + V = (26) 



The system (I20I1 - H23I1 has to be supplemented by the 
initial conditions for <fi l , Pi, a, and H. We always set 
a(t = 0) = 1, and use the first integral of motion to 

determine H{t = 0) = yj P^'+ v W+WU l ^ 

The energy density p D and the pressure p D for scalar 
dark energy are given by 



Pd 



Ekr 



V = \ Gij 



— Ekin — V — -Gi 



V&) , (27) 
V(4>) . (28) 



The total energy also includes the energy of matter, 



p T = Pm 



E, 



kin 



V 



(29) 



The f2 parameter for dark energy (matter) is given by 
the ratio of the dark energy (matter) to the total energy, 



O - Fd 

Pt 



O — Hm 
"M — 

Pt 



(30) 



Another important characteristic of the dark energy is its 
pressure-to-energy ratio defining the dark energy equa- 
tion of state p D = w D p D , 



Pl 
Pi 



Eki 



V 



E. 



kin 



V 



(31) 



The file FRWJJ2 . nb shows how we solve the Friedmann 
equations i|20[l - i|23|) f° r the N — 2 supergravity model 0] 
as discussed in the context of dark energy in [jj . Solutions 
of these equations with various initial conditions for the 
scalar fields, have an attractor behavior: the scalars even- 
tually reach their attractor values defined by a minimum 
of the potential, and the universe asymptotically becomes 
de Sitter. Note that this takes place at ifi na i w 3. How- 
ever, at this value of tfmai, when all scalars are at their 
critical values, w — — 1 and is greater than today's 
observed value Qd ~ 0.72. Nevertheless, one may still 
use this model to describe the dark energy of the present 
universe, assuming that the scalars have not reached the 
critical point yet (they are still evolving towards it). As 
we see from the plots, at ifi na i < 0.8 we have flo ~ 0.72, 
and at this time the equation of state function w(t) takes 
various values, depending on the initial conditions, which 
are mostly different from —1 but not far from it. 



V 




160 



FIG. 1: Potential of the KKLF model depending on the vol- 
ume a and the axion a 



B. A KKLT-based toy model of dark energy 

The KKLT model 4j is given by: 

With p = a + ia the kinetic term for scalars and the 
potential V(cr, a) are: 

JL[(^) 2 + (9a) 2 ] 



V = 



aAe- 2aa (A(3 + aa) + 3e aCT W cos[aa]) D 
6a 2 (T 3 



We study this model in FRW_DarkE.nb. The model has 
a number of interesting features. In some cases, when 
one starts close to the ridge, the volume first tends to in- 
crease, however, as soon as the axion starts moving, the 
volume comes back and stops there. For a long time it 
does not change, whereas the axion tends to move quickly 
towards its minimum and oscillates around it. In prac- 
tically all cases, Qd grows during the axion motion and 
reaches the value 0.72 at t » 1. During the same pe- 
riod of time w(t) remains close to -1, however it deviates 
from this value and takes various shapes, depending on 
initial conditions. We also plot the acceleration paxaxne- 

One has to keep in mind, that any realistic dark energy 
model has to accommodate today's value of the dark en- 
ergy, which means that the minimum of the potential has 
to be close to 10~ 120 . 

It is also interesting to study the late time behavior 
of the system. In FRW_LateDarkE.nb, we have solved the 
same equations with the same initial conditions up to the 
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FIG. 6: Equation of state function w(t) 



time t « 30. In most cases the system ends up at the dS 
minimum at a ~ 116 and a = 0, with w(t) = — 1. In 
other cases, both the volume and the axion have a run- 
away behavior, and the interna!7 space de-compactifies. 



IV. DISCUSSION 



FIG. 3: Volume as a function of time 




FIG. 4: Axion as a function of time 



The SuperCosmology Mathematica package presented 
in this paper proved to be useful in numerous applica- 
tions, both in our previous work as well as in the new 
models described in this paper. The interesting features 
of the new models are due to the special choice of poten- 
tials allowed in supersymmetric theories, and the non- 
canonical geometric kinetic terms. 

One of the new models studied here is the KKLT-based 
model of dark energy. The kinetic term of the model 
has an SL(2, R) -symmetry typical for string theory and 
supergravity (see [18| for earlier studies of cosmology with 
SL(2, i?)-symmetry) . One can find a change of variables 
which will bring one of the scalars to the canonical form, 

p = a + ia, a = eV^/^, however, the other one cannot 
be canonical: 



dpdp 

''{p + W 



The KKLT non-perturbative potential depends on the 



6 



volume modulus a and on the axion a and has a com- 
plicated profile. We have shown the contour plot of this 
potential as well as some trajectories of scalar fields in 
Fig. 2. 

Consider, for example, the evolution of the model from 
the point a = 130, a = 15. This initial point corresponds 
to the top left (orange) dot on the Fig. [21 The solution 
obtained numerically shows that the volume modulus a 
after some initial increase stops and waits, while the ax- 
ion evolves towards a — and oscillates around it. Then 
a moves back and eventually both fields get trapped at 
the minimum of the potential. The unusual behavior of 
the fields is explained by an interplay between the po- 
tential and non-canonical kinetic terms. From a more 
general perspective: in the second order equation for the 
scalars (|14[1 there is an extra term T l ^ k (jP(j) k in addition to 
the standard friction due to the Hubble parameter and 
also the contribution of the potential depends on the met- 
ric, G 13 That makes it hard to guess, before a numer- 
ical solution of equations is found, why some trajectories 
end up at the minimum of the potential whereas some 
other trajectories lead to the volume de-compactification. 



Our package "SuperCosmology" is intended to aid in 
the study of models related to string theory as shown in 
our examples. We hope it will prove to be useful for fur- 
ther investigations of the interface between string theory, 
supergravity and cosmology. 
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